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ABSTRACT 

We combine the estimated metallicities [Fe/H], abundances [ar/Fe], positions and motions of 
a sample of 27,500 local (7<ft/kpc<9, 0.5<|z|/kpc<2.5) SDSS/SEGUE G-type dwarf stars to 
investigate the chemo-orbital properties of the Milky Way's disk around the Sun. When we 
derive the orbital properties reflecting angular momentum, circularity, and thickness as func- 
tion of [ar/Fe] vs. [Fe/H], we find that there is a smooth variation with [ar/Fe], a proxy for age. 
At the same time, the orbital properties of the old stars with [a/Fe]>0.25 do show a transi- 
tion with [Fe/H]: below [Fe/H] ^-0.6 the orbital angular momentum decreases, and the orbits 
become significantly non-circular and thicker. Radial migration of stars into the Solar neigh- 
borhood would naturally result in a smooth variation in the orbital properties, but the latter 
old metal-poor stars form a clear challenge, in particular because a basic feature of radial mi- 
gration is that stars remain on near-circular orbits. When we next select stars on near-circular 
orbits, we indeed find besides the a-young 'thin-disk' stars a significant contribution to the a- 
old 'thick-disk' metal-rich stars. However, the remaining a-old 'thick-disk' stars on eccentric 
orbits, including nearly all old metal-poor stars, are difficult to explain with radial migration 
alone, but might have formed through early-on gas-rich mergers. We thus find chemo-orbital 
evidence that the thicker component of the Milky Way disk is not distinct from the thin com- 
ponent as expected from smooth internal evolution through radial migration, except for the old 
metal-poor stars with different orbital properties which could be part of a distinct thick-disk 
component formed through an external mechanism. 

Key words: Galaxy: abundances - Galaxy disc - Galaxy: formation - Galaxy: kinematics 
and dynamics 



1 INTRODUCTION However, a fourth, alternative model, is the in-situ forma- 

tion of the thick disk through radial migration of stars as a con- 
sequence of co-rotat i on res onance with trans ient spiral structures 
dSellwood & BinnevI [2002), bar struct ures ( Minchev & Famaevl 
2010), or orbiting satellites dOuillen et al]|200^) . This seemingly 
inevitable internal mechanism not only complicates constraining 
the contribution, if any, of external mechanisms, but as it is effec- 
tive ly a diffusion process also dilutes th e fossil record of disk galax- 
ies (freeman & Bland- Hawfhorr] |2002|) . Moreover, such a quies- 
cent internal dynamical evolution is expected to yield smooth 
changes in the properties of the disk and not distinctive thin and 
thick disk components. Still enforcing such a distinction then may 
lead to severe biases. For example, commonly adopted double- 
exponential decomposition of the light distribution yield that the ra- 
dial scale length of the thick disk component is at best weakly con- 
strained to be larger than that of the thin disk component, whereas 
a separation based on metallicity [Fe/H] and abundance [ff/Fe], 
yields that the radial scale length of t he older stars is de creasing in 
line with an inside-out disk growth dBovv et al]|2011al) . Adopting 
the same chemical separation, the resulting mass-weighted scale- 
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A thick disk, typically defined as the excess of the vertical luminos- 
ity or stellar-density profile over a sim ple exponential thin disk, is 
a gen eric feature of disk galaxies (e.g. Tsikoudi 1979; Burstein79 
1979)), including the Milky Way (e.g. lYoshiill982l ; lGilmore & Reiol 
1983). The stars in the thick disk compared to those in the thin 
disk are typically older, kinem atically hotter, more metal-poor 
and enhanced in a elements ( e.g.lFtinr mann 1998; Chiba & Beer j 
l200d : IProchaska etai] l2000l : iBensbv et al.N2005t) . Various mod- 
els for formation of a thick disk have been proposed: accretion 
of sta rs from disrupted satellites (e.g. iStatierl [l988: Aba di et ail 
2003), heating of the pre-existing thin disk by minor mergers (e.g. 
Ouinn. Hernquist & Fullagarll993l ; IVillalobos & Helmil2008h , and 
in-si tu triggered star formation during and a fter a gas-rich merger 
(e.g. IJones & Wvsdll983l : iBrook et al]|2004l) . In the paradigm of 
hierarchical structure formation all three of these external mecha- 
nisms may happen and constraining their (un)importance to thick- 
disk formation thus helps in our understanding of galaxy formation. 
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height distribution is smoothly varying an d hence opposes a dis- 
tinction in thin and thick disk components bow etaljboilbh . 

The latter authors created sub-samples of stars independent of 
the disk structure using [a/Fe] and [Fe/H] values as tags for prop- 
erties of stars that remain constant even if their orbits are changing. 
Here, we study the disk orbits in [o/Fe] vs. [Fe/H] to investigate 
whether they are also smoothly varying like the disk structure, or 
if a distinctive thick-disk component is still present, and how this 
could be understood through either of the four above thick-disk for- 
mation mechanisms. We use the SDSS/SEGUE G-dwarf sample, 
which selecte d through simple color cuts provides minimum se- 
lection biases dYannv et alJl2009l : lBovv et alj|2011al) . and includes 
reliable measurements of all six position and veloci ty com ponents 
as well as estimates of [Fe/H] and [or/Fe] i lLee et alj|2011<j) for the 
considered stars within a 1 kpc radius cylinder around the Sun. 

F or the same sample of G-type dwarf stars, IDierickx et al.l 
derive the orbital eccentricity distribution and compare 
it with that pre dicted by each of the four above thick-disk for- 
mation models dSalesetalJ 120091) . Though no firm conclusions 
are reached, the comparison seems to prefer the merger scenario, 
or radia l migration combined with another mechanism. ILee et all 
d201 lbh split the G dwarfs in a thin and thick disk subsample based 
on their bi-modal number density distribution in [a/Fe] and study 
the average chemo-kinematic properties as well as orbital ellipticity 
distribution of both subsamples to conclude that radial migration is 
important for the the thin-disk stars, but less so for the thick disk 
stars. However, it might well be the bi-modality induces a mislead- 
ing thin and thick disk separation as it is likely the natural conse- 
quence of e fficient injection of metals a t the on-set of supernovae 
type la (e.g. ISchonrich & Binnev 2009), combined with potential 
biases in the number density due to incompleteness in the SEGUE 
spectroscopic survey. 



2 SDSS/SEGUE G-TYPE DWARF STARS 

The Sloan Exten sion for Ga lactic Understanding and Exploration 
(SEGUE; lYannv et alj|2009h a s a s ub-survey of the Sloan Digital 
Sky Survey (SPSS: I York et alJ koOO) has obtained low-resolution 
(R as 2000) spectra in 3°-diameter pencil beams covering a wide 
range of sky with S > -20°. In each line-of-sight, two exposures of 
each 640 fibers per plate are made, one for bright objects and one 
for faint objects, resulting in a total of ~240,000 objects observed 
down to 20.3 mag in the g-band. Of the wide variety of type of stars 
covered, we focus on the G-type dwarf stars as they are abundant 
and have been targeted for spectroscopy with minimum selection 
biases. 

2.1 Sample selection 

We use the same selection criteria as in IDierickx et al.l d2010h to 
query the SEGUE G-dwarf sample: (i) absorption-corrected and 
dereddened color index (g - r) between 0.48 and 0.55, equal 
to the SEGUE targeting condition for G-dwarfs, (ii) color cuts 
0.6 < (;( - g)o < 2.0 and -0.1 < (r — i) < 0.4 to ensure nor- 
mal s tars, (iii) logg > 3.75 to eliminate the giant stars dYannv et al.l 
2009), (iv) E(B — V) < 0.3 to minimize effects due to uncertainty 
in extinction, (v) signal-to-noise S/N>15, and (vi) availability of 
heliocentric velocities and proper motions. Those stars that were 
not targeted as G stars but still fall within the above range of colors 
are eliminated to avoid complicated selection biases. This leads to 
a total of 27,505 G dwarfs as a starting point of our investigation. 

The correction for the incompleteness in the spectroscopic sur- 
vey, as described in the next subsection, requires sufficient stars in 
each line-of-sight so that we only consider plates with at least 100 
G dwarfs. Furthermore, we concentrate on the Solar cylinder with 
stars between 7 and 9 kpc from the Galactic center and between 0.5 
and 2.5 kpc away from the mid-plane. In the end, the sample then 
consists of a total of 14,81 1 stars. 



We correct for the latter incompleteness in the number density 
through comparison with the complete SDSS photometric survey as 
described in Section[2] together with our G-dwarf sample definition 
and derivation of their chemo-orbital properties. We then present 
maps of the average orbital properties reflecting angular momen- 
tum, circularity, and thickness as function of [ff/Fe] vs. [Fe/H] in 
Section [3] Different from the bi-modality in the number density, 
we find a smooth variation in orbital properties with [o/Fe], a proxy 
for age, as expected from internal evolution through stars migrat- 
ing in radius. In case of radial migration,the stars are expected to 
maintain their low orbital eccentricity, and indeed when select stars 
on near-circular orbits, they explain the number density distribu- 
tion as well as smoothly varying orbital properties, except for the 
old metal-poor stars. In Section [4] we discuss how current radial 
migration models seem unable to explain the latter stars with dis- 
tinct orbital properties, but how they might have formed through 
an external mechanism like early-on gas-rich mergers. This leads 
us in Section [5] to the conclusion that the Milky Way's thick disk 
has smoothly evolved from the (single-exponential) thin disk prob- 
ably through radial migration, except for the old, mainly metal-poor 
stars with distinct orbital properties which could be part of a dis- 
tinct thick-disk component formed through an external mechanism. 

Throughout we adopt 8 kpc for the Sun's distance to the 
Galactic center, and 220 km s -1 for the circular velocity of the local 
standard of rest (LSR). 



2.2 Star count correction and number density 

The spectroscopically targeted stars are randomly selected from 
a color-selected photometric sample, so that their star count N spcc 
can be corrected for by comparing with the photometric star count 
Afphot in the same volume. Along a line-of-sight pencil beam, the 
star count is thus a function of colors, e.g., (u - g)o and (g - r)o, 
and distance d. Distances to individual stars are estimated accord- 
in g the photometric c olor-metallicity-absol ute-magnitude relation 
of llvezig et ail d2008l ; their eq. A7). While iKlement et al.l J2009b 
suggest that based on the analysis of 11 stellar cluster this pho- 
tometric method may systemati cally underestimate the distance by 
about 3 percent. lAn et aljd2008l) find that based on the distances de- 
rived from stellar isochrones the photometric distances are about 10 
percent larger. Fortunately, the distances uncertai nties seem to hav e 
little influence on the stellar density (see also ljuric et al.l 1 2008). 
and to be independent of chemical properties (see also lBovv et al.l 
1201 1 ah . so that they will not affect our chemical-orbital analysis. 

Assuming that within a small range of (u - g)o, (g - r)o, and d 
the spectroscopic selection is indeed random, i.e., potential depen- 
dencies on other properties are negligible, the star count correction 
of the spectroscopic sample follows from multiplying with 

w = N phot ((u - g) , (g - r) , d)/N spec ((u - g)o, (g - r) , d). (1) 

When computing average properties from the spectroscopically tar- 
geted stars, the inclusion of the latter weights ensures that stars 
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Figure 1. Equatorial map of the SDSS/SEGUE G-type dwarf stars in our sample with radius 7 < S/kpc < 9 from the Galactic center and height 0.5 < 
Izl/kpc < 2.5 away from the Galactic mid-plane. The colors with corresponding values given by the bar on the right side indicate the ratio of photometrically 
to spectroscopically observed G-dwarfs with same colors (u - g)o and (g - r)o and distance. These ratios are used as weights to correct the star counts of the 
spectroscopic SEGUE G-dwarf sample, which is incomplete, especially toward lower latitudes as indicated by the redder colors. 



which are spectroscopically under-sampled are boosted through 
w > 1, and stars that are spectroscopically over-sampled are sup- 
pressed through w < 1 . 

In Figure [JJ the stars in each of the 162 SEGUE plates are 
color coded with the value of w. The stars toward the Galactic plane 
have a relatively high weight w > 1 (reddish color) as a result of the 
high density of stars with respect to the number of fibers available 
per plate for spectroscopic follow-up. Toward the Galactic poles the 
density is such that typically spectra are taken of all stars, so that 
little correction is needed leading to weights w ~ 1 (blueish color). 
To avoid large and potential over-corrections, we select stars with 
weights 0.5 < w < 10, reducing the total number of G dwarfs to 
13,084. 

Even though our sample of stars is restricted within a narrow 
range of radii of 7 < R < 9 kpc, the variation with radius R is 
non-negligible when deriving the stellar numb er density v as func- 
tion of height \z\ above the mid-plane. Whereas lBovv et all ( l2011al) 
perform the non-trivial structural fit in R and |z| simultaneously, be- 
low we fit only for the exponential scale height while adopting an 
exponential in radius with fixed scale lengths form these authors. 
Specifically, we use scale lengths of 1.9 and 2.2 kpc for stars with 
metallicities below and above [Fe/H]=-0.6, although the results are 
insensitive to the precise values, as long as the exponential depen- 
dence with radius is taken into account. 



2.3 Metallicities, abundances and space motions 

Measurements of metallicity [Fe/H] and abundance [ff/Fe] of the 
G dwarf stars are obtained from the S EGUE Stellar Parame - 
ter Pipeline (sspp) table in SDSS DR7 jAbazaiian etail 120081) . 
Through comparisons with five stellar clusters, the typical uncer- 
tainty in [Fe/H] is estimated at 0.13 dex, but also systema tically un- 
deres timated by about 0. 1 dex for solar-abundance stars dLee et all 
2008). The latter skews the metallicity distribution function (MDF) 
around [Fe/H]=0 to lower values, but as we are mainly interested 
in the kinematics (well) below solar, this does not play a role in 
our analysis. Through comparisons with estimates from higher- 
resolution spectra, the typic al uncertainty in [ff/Fe] is estimated 
at 0.1 dex dLee et al]|2011ah . for stars with effective temperature 
4,500 < 7/rff < 7,000 (in Kelvin), surface gravity 1.5 < logg < 
5.0, and metallicity -1.4 < [Fe/H]^ +0.3, covering most of the G 



dwarfs in our sample. This 0.1 dex uncertainty is based on SEGUE 
spectra with S/N>20, but our more inclusive S/N>15 will not in- 
crease the uncertainties on [or/Fe] by much. 

The measured line-of-sight velocities, often confusingly also 
called radial velocities, and proper motions of the stars are trans- 
formed into th e three velocity components along cylindrical coordi- 
nates (see e.g. ljohnson & Soderblom| [T987; Kle ment. Fuchs & Rixl 
2008), namely radial velocity v R , azimuthal or rotational veloc- 
ity v^, and vertical velocity v z . In case of v^, positive values re- 
fer to the same direction of the rotation of LSR, and stars with 
= 220kms~' move at the same speed as the LSR while those 
with = km s have no (azimuthal) rotation. We adopt for the 
Sun's peculiar velocity relative to the LSR the common values of 
(10.00, 5.25, 7.17) kms~' in the radial, azimuthal and vertical direc- 
tion, respec tively, and ev en though the azimuthal component might 
be revised (Binney 2010) it does not affect our analysis. 



2.4 Orbital properties 

Given a star's current position and space motions, we com- 
pute its orbit in a three-co mponent Galactic potential Q(R, z), in- 
cludi ng a Hernquist b ulge jHernquist|[T99fj|) . a Miyamoto-Nagai 
disk ( Miyamoto & Nagai 1975), and a NFW dark-matter halo 
i Navarro, Frenk & White 1996), with parameters adopted from 
Gome zet alj d2010l) . We next obtain the guiding radius R g of 
each star by solving the equation d<b/dR = J./R 3 in the mid- 
plane (z = 0), with J z = Rv<fi the angular momentum paral- 
lel to the symmetry z-axis. The orbital eccentricity is defined as 
e = (R apo ~ ^pen)/(^a P o + ^peri), with R peri and R apo the peri-center 
and apo-center radius of the computed orbit, respectively. We also 
record v z max as the maximum velocity with which the stellar orbit 
is crossing the mid-plane as an indicator of the vertical energy or 
thickness of the orbit. Finally, as an alternative measure of the or- 
bital eccentricity we compute JJJ C , the ratio of the z-component of 
the orbital angular momentum to the maximum angular momentum 
J c = R c v c when the orbit is circular. The latter radius R c follows 
from solving E = O + v^/2 in the mid-plane, with E the orbital 
energy and v c the circular velocity defined as v 2 c = Rd<b/dR for 
z, = 0. 
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Figure 2. Star-count-corrected number density of G dwarf's in abundance [»/Fe] vs. metallicity [Fe/H]. Left: The number density of all G dwarfs shows a 
clear bi-modality that naturally inspires a separation around [ti'/Fe]=0.25, indicated with a horizontal dashed line, into a-young 'thin-disk' stars and a-old 
'thick-disk' stars. The vertical dashed line at [Fe/H]=-0.6 indicates roughly the transition in the orbital properties of the o--old stars shown in Figure [3] In 
this way the G dwarfs are split into the three groups indicated with the Roman numerals, while the bottom-left corner contains only relatively very few stars. 
Middle and right: The number density of G dwarfs which are on (near-)circular orbits with 7-// c ,>0.85, and on eccentric orbits with / Z /7 C ,<0.80, respectively. 
Here, /. is the angular momentum momentum parallel to the symmetry z-axis, and J c the (maximum) angular momentum of a circular orbit. As expected 
nearly all of the a-young, metal-rich stars in group I are on near-circular orbits. However, also a significant fraction of the a-old, metal-rich stars in group II 
are on near-circular orbits, while most of the ff-old, metal-poor stars in group III are on eccentric orbits. 



3 CHEMO-ORBITAL ANALYSIS 

We use the star-count-corrected spectroscopic G dwarf sample 
to study their chemo-orbital properties. Specifically, we investi- 
gate the orbital angular momentum, circularity and thickness of 
G dwarfs as function of their metallicity [Fe/H] and abundance 
[a/Fe]. We use the latter as an age indicator, and hereafter refer 
to stars with abundances [a/Fe] >0. 25 as a-old stars, and stars with 
[a/Fe]<0.25 as a-young stars. The split is initiated by a bi-modality 
in the number density distribution of G dwarfs around [a/Fe]=0.25, 
but as argued below does not necessarily imply the presence of a 
bi-modality in age. 



3.1 Smooth variation orbital properties 

The left panel of Figure [2] shows the map of the star-count cor- 
rected number density of G dwarfs in [a/ Fe] vs. [Fe/H], O verall 
the distribution is similar to that found bv lLeeetal.U2011bh . even 
though we include a star-count correction, adopt a different selec- 
tion of G dwarfs, and are using (slightly) different measurements of 
[ff/Fe] and [Fe/H]. A bi-modality is clearly visible that naturally in- 
spires a separation around [a/Fe] =0.25 (indicated with a horizontal 
dot-dashed line) into low er- [g/Fe] thin-disk stars and high- [a/Fe] 
thick-disk stars, similar to lLee et al.l feOllbh . 

However, even though [a/Fe] is an age indicator with higher 
[a/Fe] corresponding to older stars, the bi-modality in [a/Fe] does 
not necessarily imply a same bi-modality in age. Instead the bi- 
modality could be the result of the on-set of supernovae type la 
that efficiently inject metals but not a-elements into the interstel- 
lar gas, causing a sudden drop in [a/Fe] for the newly formed 
stars. Moreover, the local G dwarfs have been selected in the ra- 
dial range 7 < R < 9kpc, so that the spread in radii where the 
stars are born is at leasy 2kpc. Because in general star forma- 

1 The range in birth radii is likely larger due to the presence of stars on 
eccentric orbits and in particular in case of radially migrated stars. 



tion is higher at smaller radii due to higher gas surface density, the 
drop in [ff/Fe] due to the on-set of supemovae type la should occur 
at higher [Fe/H] at smaller radii. The apparent bi-modality in the 
number density over a range of metallicities, might thus be the nat- 
ural result of the mix of supernovae type la enrichments that hap- 
pened at a range of radii. Consequently, while [a/Fe] has changed 
abruptly, the star-formation history and subsequent evolution of the 
stars might have been smooth throughout. In this case, we expect 
their orbital properties to change also smoothly with [a/Fe] and 
[Fe/H]. 

Figure [3] shows maps of the average orbital properties of the 
G dwarfs in [a/Fe] vs. [Fe/H], On the left are maps of the median 
azimuthal rotational velocity and median guiding radius R g , both 
indicative of the orbital angular momentum. In the middle are maps 
of median eccentricity e and relative-to-circular angular momen- 
tum J z /J c , both indicative of the orbital circularity. On the right 
are maps of vertical velocity dispersion a- and average squared 
maximum vertical velocity (v 2 max ) 1/2 , both indicative of the orbital 
thickness. In all maps the average orbital properties are indeed vary- 
ing smoothly with [a/Fe]. 

The only apparent transition is with [Fe/H] for the a-old stars: 
below [Fe/H]^-0.6 (indicated with a vertical dot-dashed line) the 
orbital angular momentum decreases, and the orbits become sig- 
nificantly non-circular and thicker. Particularly striking is this tran- 
sition in (v 2 max ) 1/2 (bottom-right panel): (v 2 max ) 1/2 smoothly in- 
creases with [a/Fe] for the metal-richer stars, but below [Fe/H]=:- 
0. 6 there is a rather abrupt increase to < v? max > 1 12 a 1 20 km s~ 1 , seem- 
ingly independent of [a/Fe] (>0.25). 

This transition in [Fe/H] for the a-old stars could be the result 
of probing increasing heights \z\ above the equatorial plane with de- 
creasing [Fe/H]. In other words, if instead of including stars at all 
heights, we only consider stars within a small range of \z\, would 
the transition with [Fe/H] disappear? To maintain a sufficient num- 
ber of stars per bin in we adopt the splits at [a/Fe] =0.25 and 
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Figure 3. Average orbital properties of G dwarfs in [a/Fe] vs. [Fe/H]. As in Figure[2] the stars are split by a horizontal dashed line at [a/Fe]=0.25 and a vertical 
dashed line at [Fe/H]=-0.6 in three groups indicated by Roman numerals. Left column: Maps of median azimuthal rotational velocity Vd, and median guiding 
radius R g , both indicative of the orbital angular momentum. The anti-correlation with [Fe/H] for the ti-young, metal-rich stars (group I) disappears at higher 
[a/Fe], but the the inversion toward a positive correlation is actually only evident for the tr-old, metal-poor stars (group III), while the a-old, metal-rich stars 
(group II) show no apparent correlation. Middle column: Maps of median eccentricity e and relative-to-circular angular momentum (J-/J c ), both indicative of 
the orbital circularity. Whereas at higher metallicity, stars with both lower and higher [a/Fe] (group I and group II) are on close-to-circular circular orbits, only 
the metal-poor stars (group III) are on significant non-circular orbits. Right column: Maps of vertical velocity dispersion cr. and average squared maximum 
vertical velocity (v?™^) 1 ' 2 , both indicative of the orbital thickness. While for metal-rich stars the orbital thickness smoothly increases independent of [Fe/H], 
the orbits of the metal-poor stars all seem to be significantly thicker. 



[Fe/H]=-0.6 to arrive at three groups^ indicated by the Roman nu- 
merals in Figures [2] and [3] group I of a-young, metal-rich stars, 
group II of a-old, metal-rich stars, and group III of ff-old, metal- 
poor stars. Figure [4] shows the same average orbital properties but 
now as function of height for these three groups. The offset in 
average orbital properties between group I (magenta curves) and 
group II (green curves) are expected from the difference in [o'/Fe] 
and thus age. However, group II (green curves) and group III (or- 
ange curves) are both a-old, but still show a clear offset, which 
moreover remains nearly constant with height. This means if we 
would select a narrow range in z the transition for the a-old stars 
with [Fe/H] should remain. In other words, the marginalization over 
height has not affected the overall trends seen in the maps in Fig- 
ures[2]and[3] including the smooth change with [ff/Fe]. 



3.2 Radial migration 

This smooth variation of the orbital properties with [a/Fe] and 
hence age seems consistent with recent numerical implementations 



ot m-situ tormation ot the thick disk through radial migration (e.g. 
Schon rich & Binn ev 2009: Eoebman et alj 201 lh . However, these 



A fourth group of a-young, metal-poor stars contains too few stars. 



radial migration models have difficulty explaining various proper- 
ties of the o--old G dwarfs . 

First, lLoebman et al] J201 lh find that the number density dis- 
tribution of their simulated Solar Neighborhood is not bi-modal in 
[a/Fe] (see their Fig. 12) opposite to the clear bi-modality for G 
dwarfs (see left panel of Figuref2](. Still, this might partially be the 
result of the adopted chemical model in their simulation, as SB 10 
on the other hand do seem to find a depression in the number den- 
sity distribution around intermediate abundance values (see their 
Fig. 3). 

Second, Sc honrich & Binnevl fc009h find that the mean rota- 
tion velocity (v^) is anti-correlated with mean metallicity [Fe/H] 
similarly for all their simulated stars (see their Fig. 4), while 
the anti-correlation clearly weakens and even becomes inverted 
for the a-old G dwarfs (see top-left panel of Figure O. How- 
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Figure 4. Average orbital properties of G dwarfs as in Figure[3]but now as function of height \z\ away from the equatorial plane. The three different curves in 
each panel correspond to the three groups in [a/Fe] vs. [Fe/H] as indicated in Figure|2] connected magenta circles for cr-young, metal-rich stars in group I, 
connected green diamonds for cr-old, metal-rich stars in group II, and connected orange triangles for ff-old, metal-poor stars in group III. In all panels, all 
three groups show a close to linear trend with similar slopes but different intercepts, independent of height. While such offsets might be expected between the 
a-young stars (magenta) and a-old stars (green and orange) because of difference in age, the systematic offset between the cr-old, metal-rich stars (green) and 
ff-old, metal-poor stars (orange) at the same [a/Fe] imply another explanation than age. Of the orbital properties, <x- as second velocity moment is the most 
sensitive to outliers and number of stars per bin in |z|. While we take into account both effects, fluctuations with height remain present, but the offsets between 
the three groups are still evident. 



ever, Schonrich & Binney ( 2009) select stars from their simulation 
within 100 pc from Sun while the sample of G dwarfs only starts be- 
yond ~500 pc, at which point (in height above the equatorial plane) 
the or-old stars with diluted rotation-metallicity anti-correlation 
signifi cantly increase in relative numbers. Indeed, lLoebman et al.l 
( 120 111) probing larger height find that the anti-correlation weakens 
with age and hence [ff/Fe] (see their Fig. 10). 

Third, the ff-old, metal-poor stars (in group III) seem to form 
the biggest chal lenge for radial mig r ation mode ls. In the sim- 
ulation s of both Schonrich & Binnevl d2009h and lLoebman et al.l 
d20 1 lh most of the stars that radially migrated into the Solar Neigh- 
borhood were born in the inner regions, and hence, if any, will 
be more metal-rich than the stars born in-situ around the same 
time. Perhaps even more problematic is that on average the a- 
old, metal-poor stars are on quite eccentric orbits (see middle pan- 
els of Figure [3}, whereas a basic property of radial migration is 
that stars remain on near-c ircular orbits when they move in radius 
dSellwood & Binnevll2002l) . 

As we have computed the orbital circularity for each G dwarf 
star, the latter naturally leads to question what happens if we select 
only the stars on near-circular orbits. Instead of eccentricity e, we 
use the relative- to-circular angular momentum J z /J c , which is a 
dynamical measure of orbital circularity that is less sensitive to the 
adopted underlying mass model for the Milky Way. 

Even if at their birth radii the stars move on circular orbits 
with J z /J c =l, the current distribution in J z /J c at the Solar Neigh- 
borhood is expected to extend toward values below unity, both be- 
cause of stars scattering to more eccentric orbits, as well as due to 
measurement errors, mainly in their distance and proper motions. 
Indeed, selecting the a-young, metal-rich stars of group I which are 
expected to be on near-circular orbits, yields a distribution in J z /J c 
(blue curve in bottom-left panel of Figure[9} that is strongly peaked 
toward unity but also has a significant tail to lower values. In what 
follows, for simplicity we select stars with J Z /J C >0.S5 to be a on 



near-circular orbits. Below in Section |4~3l we argue that our results 
are robust against varying the latter lower limit on J z /J c , or when 
we use the J z /J c distribution of group I to assign per star a prob- 
ability that it is on a similar near-circular orbit as these 'thin-disk' 
stars. 

3.3 Stars on near-circular orbits 

The middle panel of Figure [2] shows the number density distribu- 
tion in [or/Fe] vs. [Fe/H] for G dwarfs on (near-)circular orbits with 
7 Z /7 C >0.85. For the young, metal-rich stars (group I), the distribu- 
tion is nearly unchanged from the number density of all G dwarfs 
(left panel), as expected for thin-disk stars that mostly move on or- 
bits that are not far from circular. At the same time, there remains 
a significant fraction of old, metal-rich stars (group II), but only 
very few old, metal-poor stars (group III). The number of old stars 
has decreased relative to the younger stars such that there is not 
anymore a clear bi-modality in [ff/Fe]. 

The [a/FeHFe/H] maps in Figure [5] reveal that there are also 
significant changes in the average orbital properties of the a-old 
stars when only those on near-circular orbits are selected: the aver- 
age orbital angular momentum (left panels) has significantly in- 
creased, while the average orbital thickness (middle panels) has 
significantly decreased. Furthermore, the transition with [Fe/H] has 
disappeared, with now orbital angular momentum and thickness on 
average the same for all [Fe/H] values. 

The top-right panel of Figure [5] shows the number density as 
function of height |z| away from the equatorial plan^l for the same 
G dwarfs on near-circular orbits. The decrease in slope at larger 
height is well fitted by a double exponential, with corresponding 

3 The range in radii 7 < R < 9 kpc has been taken into account through a 
exponential correction term in radius as described in Section [Z2l 
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Figure 5. Average properties of G dwarfs which are on (near-)circular orbits with J-/ J c ,>0.85. As in Figure[3] the average orbital properties in the first and 
second column are indicative of the angular momentum and thickness of the orbits in [a/Fe] vs. [Fe/H]. While va and R g are (weakly) anti-correlated with 
[Fe/H] only for the ff-young stars, cr. and (vf™*) 1 ' seem to be independent of [Fe/H] and to increase smoothly with increasing [a/Fe]. The third column 
shows the number density v and vertical velocity dispersion cr z as function of height |z| away from the equatorial plane. The black solid curve shows the best-fit 
double exponential fit to the number density profiles, with corresponding best-fit scale heights and relative fraction of the second exponential at z = indicated. 
The vertical dispersion profile is nearly constant (isothermal) with average value indicated. 



scale heights of zo = 281 pc and zi = 901 pc in full agreement 
with e arlier determina tions of the thin and thick disk scale heights 
(e.g. Ijuric et alj [2008). The vertical velocity dispersion tr. in the 
bottom-right panel varies only mildly around the average value of 
about 27kms -1 , indicating that the G dwarf stars on near-circular 
orbits are dynamically close to a simple isothermal model in height. 

3.4 Stars on eccentric orbits 

To have a look at the average properties of the G dwarf stars that 
are on eccentric orbits, we select all stars with J z /J c <0.S0. As for 
selecting the stars on near-circular orbits this upper limit in J z /J c 
is chosen for simplicity, but we argue below in Section [4731 that 
our results do not depend on the details of this selection. The num- 
ber density distribution in [a/Fe] vs. [Fe/H] in the right panel of 
Figure [2] shows that, as expected, nearly no G dwarfs on eccentric 
orbits are cr-young, metal-rich stars (group I), that they contribute 
relatively little to the ff-old, metal-rich stars (group II), but that they 
dominate the o--old, metal-poor stars (group III). 

The [a/Fe]-[Fe/H] maps in Figure [6] reveal that the average 
orbital properties show no significant trends with neither [a/Fe], 
nor with [Fe/H]. The average orbital angular momentum and thick- 
ness are respectively much lower and higher than G dwarfs on 
near-circular orbits. For example, the average azimuthal rotation 



velocity (v^) is about half, and the vertical velocity dispersion <x- 
about twice that of the o--old stars on near-circular orbits. The top- 
right panel shows that the vertical number density is well fitted by 
a single exponential, with scale height zo = 865 pc that is fully 
consistent with the larger scale height of the double exponential 
fit to stars on near-circular orbits and thus als o with earlier de - 
terminations of the thick-disk scale height (e.g. Ijuric et ali 2008). 
The local eccentric-to-near-c ircular dens i ty nor malization is about 
3.5%, which is smaller than Ijuric et al. (2008) who derived 12% 
for thick disk. It implies that the thick disk defined by the double- 
exponential vertical profile must contain the stars both on eccentric 
and near-circular orbits, which are probably originated from dif- 
ferent mechanisms. The vertical velocity dispersion in the bottom- 
right panel again varies only mildly around the average value of 
about 58kms~', indicating that the G dwarfs on eccentric orbits 
are dynamically also close to a simple isothermal model in height. 



4 DISCUSSION 

When dividing the G dwarf stars based on the circularity of their 
orbits, we found that the G dwarfs on near-circular orbits capture, 
as expected, the a-young, metal-rich 'thin-disk' stars. However, G 
dwarfs on near-circular orbits also make up a significant fraction of 
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Figure 6. Same as Figure[5] but for G dwarfs which are on eccentric orbits with J z /J c ,<0.H0. Stars on eccentric orbits are mainly a-old, metal-poor stars (see 
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as expected, the exponential scale height of these a-old 'thick-disk' stars is larger than that of the a-y oung 'thin-disk' stars in Figure[5] the significantly lower 
average guiding radius would implies a smaller scale length for the older stars (cf. lBovv et alj201 lal) , in line with an inside-out disk growth. 



the a-old, 'thick-disk' stars, and explain the average orbital prop- 
erties of the metal-rich stars that vary smoothly with [a/Fe]. The 
remainder of the a-old, mainly metal-poor stars are covered by G 
dwarfs on eccentric orbits, with significantly different average or- 
bital properties which in turn explains the apparent transition with 
[Fe/H] when the a-old stars are taken together. Does this imply that 
the thick disk is made of a-old stars that have (at least) two different 
origins, reflected in their orbital properties? 



4.1 Origin of a-old stars on near-circular orbits 

It seems very plausible that radial migration is able to explain 
the properties of the G d warfs on near-circular orbits. Since 
ISchonrich & Binnevl |2009» analyze their radial migration simula- 
tion for stars within 100 pc from the Sun, while the sample of G 
dwarfs only star ts beyond ~500pc, w e focus here on the simula- 
tion results from lLoebman et all J201 lh . 

First, the absen ce of a bi-modality in [a/Fe] found by 
lLoebman et alj J201 it their Fig. 12) is not anymore inconsistent 
with the number density of G dwarfs after only considering those 
on near-circular orbits (middle panel of Fig. f2J- Second, the re- 
moval of most of the a-old, metal-poor stars also elevates the 
big challenge to have to explain those metal-poor stars with ra- 
dial migration that predominantly moves stars out from a relatively 
metal-richer inner region. Third, whereas for all G dwarfs together 
the anti-correlation between average azimuthal rotational velocity 



(v^) and metallicity [Fe/H] even got inverted at high [a/Fe] values 
(top-left panel of Figure [3}, considering only G dwarfs on near- 
circular orbits the anti-correlation weaken s for increasing [q/Fe] 
(top-left panel of Figure [5), consistent with lLoebman et all J20 1 ll ; 
their Fig. 10) 

The weakening of the rotation-metallicity anti-correlation 
with age is further illustrated in the left column of Figure [7] Per 
bin in [a/Fe] values, increasing from top to bottom panels as indi- 
cated, the rotation versus metallicity for G dwarfs with J-/J c >0.&5 
are plotted. In each panel the correlation coefficient of the indi- 
vidual measurements (dots) as well as the slope of the linear fit to 
the median values per bin in [Fe/H] (open circles) are given. The 
sign and amplitude of both the correlation coefficient and the fitted 
slope show how the rotation-metallicity anti-correlation weakens 
at higher [a/Fe], and even disappears fo r the older stars, consistent 
with the radial migration simulations bv lLoebman et al. 

The middle column of Figure [7j shows the vertical number 
density fitted by a single exponential, with resulting scale heights 
that are smoothly increasing with [a/Fe]. The vertical velocity dis- 
persions in the right column are only mildly varying around average 
values that are also smoothly increasing with [a/Fe], W hile this sur- 
prisingl y simple exponential behavior is investigated in lBovv et al.1 
u20 1 1 af) and the near-isothermalilty will be studied in more detail in 
forthcoming papers, it will be interesting to see if radial migration 
models predict similar vertical exponential density and flat veloc- 
ity dispersion profiles. At the same time, a more direct comparison 
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with the [(i/Fe]-[Fe/H] maps of average orbital properties shown 
in Figure [5] are necessary to (dis)prove that radial migration is the 
origin of (most of the) cr-old stars on near-circular orbits. 

4.2 Origin of cr-old stars on eccentric orbits 

While radial migration might be the main formation channel for the 
a-old stars on near-circular orbits, we already argued (Section[3]2j 
that radial migration most likely cannot be the origin of the a-old 



stars on eccentric orbits. As mentioned in SectionQ] alternative for- 
mation channels are (i) accretion of stars from disrupted satellites, 
(ii) heating of an existing thin disk by minor mergers, and (iii) in 
situ formation during/after a gas-rich merger. 

The fact that G dwarfs on eccentric orbits are mainly ff-old, 
metal-poor stars (right panel of Figure^ suggests that their forma- 
tion happened at high redshift and before the on-set of supernovae 
type la. Still, in the case of accretion from dwarf satellites, the ob- 
served dwarf galaxies with such low metallicities are also those 
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with the lowest stellar mass jDekel & Wooll2003l) . As a result, the 
dwarf satellites are expected to have had a slow star formation rate 
with relatively few supernovae type II, which in turn mak es the a- 
(pre-)enhancement difficult to explain (Ruchti et al. 2 0101) . 

In the scenario of heating by more massive satellites, the 
heated stars might have the appropriate chemical properties, but 
their orbital properties are incompatible with many of the a-old 
stars on eccentric orbits. For example, the top-left panel of Figure[8] 
shows that the azimuthal rotational velocity for most of the stars 
on eccentric orbits (red histogram) is more than 100 km s~' lagging 
the LSR (at 220km s -1 ), while such lar ge lags are neve r reached in 
simulations of satellite heating (see e.g. u et alj|201 lh . A way out 
seems to be cumulative heating through multiple minor mergers, 
but then we expect to observe some satellites on low eccentric or- 
bits around the Milky Way which are not (yet) disrupted at present 
days, and we should observe the debris from the disrupted satel- 
lites on high eccentric orbits in the disk. Neither such kind of dwarf 
galaxies on low eccentric orbits have be en found by now, nor the 
expected bump of high eccentric orbits dSales et alj|2009f) is seen 
in the eccentricity distribution in the top-middle panel of Figure[8] 

In case of a gas-rich merger at high redshift, it is expected 
that the disk from which the stars form is not s mooth and thin, bu t 
clumpy with already a significant dispersion dBrook et alj|2004l) . 
Together with further blurring over the long time since the merger, 
it might well be that if any correlation with [cr/Fe] or [Fe/H] existed 
it is washed out by now. Even more so, an early-on merger origin 



from an already hotter disk might help explain the apparent jump 
between the average orbital properties of G dwarfs on near-circular 
and eccentric orbits that can be seen by comparing the maps in 
Figures [5] and [6] For example, in case of the a--old stars around 
[Fe/H]=-0.6, going from the stars on near-circular orbits to those on 
eccentric orbits, the median v<j halves from ~200 to ~100 kms -1 , 
while cr z doubles from ~30 to ~60 km s -1 . 

In case of both accretion and heating not only satellite stars 
are expected to be accreted but also the dark matter in which the 
satellites were embedded. The dark matter will settle into a disk, 
which we might be able to detect through its gravitation al effects 
on the surrounding baryonic matter (e.g. Read et alj2008T) . As long 
as such a so-called 'dark disk' remains undetected, we believe the 
chemo-orbital evidence is in favor of a merger origin for the a-old, 
thick-disk stars on eccentric orbits. 

Using a similar sample of G dwarfs, but concentrating on 
the eccentricity d i stribution as a thick-d isk formation diagnostic 
dSales et alj|2009l) . iDierickx et al.l d2010|) also found that both ac- 
cretion and heating with too many predicted accreted stars on high- 
eccentric orbits are less likely than a merger origin. At the same 
time, they found that radial migration alone cannot explain the ec- 
centricity distribution of all G dwarfs together, fully in line with our 
findings of a mixed origin for the Milky Way's thick disk. 

Note that contamination of the a-old stars on e ccentric or- 
bits by halo stars is very small (see also section 5.2 of iBovv et al.l 
1201 la ), and does not influence our analysis. 
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Figure 9. Alternative selection of stars on (near-)circular and eccentric orbits, and the resulting star-count corrected number density of G dwarfs in [ff/Fe] vs. 
[Fe/H]. Top and middle row; Similar hard-cut in orbital circularity J z /J c as in Figure[3] but with different limits as indicated at the top of the two rows. Bottom 
row: The left panel shows in black the distribution of J z /J c for all stars in black, in blue that of the a-young, metal-rich stars of group I, re-scaled so that the 
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yields per star a probability that it is on a near-circular orbit given its value for J z /J c , resulting in the map for circular orbits in the middle panel, while the map 
for eccentric orbits in the right panel is the difference with respect to the map of all orbits. 
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4.3 Alternative chemo-orbital selection? 

Instead of eccentricity, we used the relative-to-circular angular mo- 
mentum J.J] C as a dynamical measure of orbital circularity that is 
less sensitive to the adopted underlying mass model for the Milky 
Way. However, from Figure [8] (top-middle panel) it is clear that 
selecting on eccentricity will yield similar results. We also tried 
different lower and upper limits of J,/J c for the selection of near- 
circular and eccentric orbits, but found that when we stayed within 
the reasonable range O.S<J./ J c <0.9 (top and middle rows in Fig- 
ure^, there are no significant changes in the results. 

Alternatively, instead of a hard cut in J z /J c , we also investi- 
gated assigning to each star a probability that it is on a near-circular 
orbit given its value for J z /J c . To this end, as shown in the bottom- 
left panel of Figure(8] we assign part of the distribution in J z /J c for 
all stars (black curve) to near-circular orbits (blue curve) and the 
remainder to eccentric orbits (red curve). The distribution in J z /J c 
for near-circular orbits is simply a re-scaling of the corresponding 
distribution in J z /J c for the a-young, metal-rich 'thin-disk' stars 
(group I). Adopting this probability per star, yields maps of num- 
ber density in [a/Fe] vs. [Fe/H] for stars on circular and eccentric 
orbits (bottom-middle and bottom-right panels of Figure[8](, which 
are similar to those based on hard cuts in J z /J c . However, adopt- 
ing these probabilities will somewhat blur the clear distinction in 
average orbital properties seen in Figures [5] and [6] for an adopted 
hard cut, because stars on near-circular orbits also contribute to the 
properties of the eccentric orbits through the tail in the distribu- 
tion of J-JJ C toward low values, and vice versa. Instead, one should 
use the probabilities to first classify each star into groups, and only 
then derive average orbital within each group. Even though this 
will improve the hard-cut classification adopted in this paper for 
simplicity, we believe our results on the mixed origin of the a-old 
'thick-disk' stars to be robust. 

Selection on orbital circularity is a natural consequence of the 
basic property of radial migration that stars remain on nearly cir- 
cular orbits when they move in radius. It implies that the angular 
momentum has been changed during the radial mixing, making 
or Rg overlap with those of pre-existing old metal-poor stars. In- 
deed, the angular momentum distribution for stars on circular or- 
bits overlaps with those on eccentric orbits as shown in left panels 
of Figure [8] For this reason, it would be not feasible to select stars 
based on angular momentum related parameters. On the other hand, 
we see from the bi-modal distributions in the right panels of Fig- 
ure[8] that the orbital thickness might be a good alternative to split 
the stars with different origins; o~ z is not practical as it is based on 
an ensemble of stars, whereas v ;jTlax can be inferred for each indi- 
vidual stars. Since the bi-modality of the distribution of v z raax is in 
line with the split in J z /J c , it is expected that it will yield nearly 
identical results in orbital properties. 



5 CONCLUSIONS 

Since iGilmore & Reid Jl983h fitted the vertical stellar density pro- 
file in the Solar neighborhood with a double-exponential, the Milky 
Way's stellar disk is thought to consist of a thin and thick compo- 
nent with distinct spatial as well as kinematical and chemical dis- 
tributions. However, separating the two populations in kinematics 
or chemical properties of the stars is difficult if not ambiguous. 

For instance, when stars in the Solar Neighborhood are se- 
lected to be thick-disk stars if their velocities are significantly away 
from the local standard of rest, also stars will be included that orig- 
inated from the inner disk and reached the Solar Neighborhood 



because of their eccentric orbit and/or via radial migration. This 
means that the kinematically selected thick-disk is 'contaminated' 
with stars born at smaller radii with higher metallicities, some even 
more metal-rich than the kinematically selected thin-disk stars that 
mostly originate from around the Solar radius. Moreover, a large 
fraction of these contaminating stars from the inner disk will have 
formed after supernovae type la enrichment and hence not only 
make the kinematically selected thick-disk to extend to metallic- 
ities beyond solar but at the same time also to lower a-abundances, 
i.e., introducing a 'knee' in the abundance-metallicity distribution. 
In addition, the latter 'knee' will be shifted to higher metallicities 
so that the abundance trends for the kinematically selected thin and 
thick disk are separated. These effects are indeed seen in studies of 
samples of local stars with detailed abundance meas urements that 
are ki nematically divided into thin and thick disk (e.g. Bens bv et al.l 
2005). Drawing conclusions on the thick-disk formation based on 
these findings seems, however, questionable as they depend on the 
kinematical selection tha t inevitably add a-young, m etal-rich stars 
to the thick disk (see also lSchonrich & Binnevll2009t) . 

Whereas kinematical properties change all the time, chemical 
composition is a permanent feature of a star and intimately con- 
nected to time and place of its birth. Even so, selection based on 
metallicity is ambiguous because stars of different birth radii and 
hence different metallicities contribute to the Solar neighborhood, 
giving rise to a board metallicity distribution, even at fixed age. 
Moreover, even if the thick disk is on average older than the thin 
disk, stars with a range of ages and hence metallicities can con- 
tribute to both, so that in the end a large overlap in metallicities 
between thin and thick disk stars is expected. Better is to use a- 
abundance which is relatively independent of birth radius and more 
directly linked to age, but at the same time harder to measure and 
still only a proxy for age. For example, the bi-modality in the num- 
ber density distribution of G dwarf stars with [a/Fe] (see also left 
panel of Figu re [H inspi r ed the chemical separation into thin and 
thick stars by iLee etal (2011b). However, such a bi-modality in 
[a/Fe] could well be the natural consequence of the supernovae 
type la enrichment physics and not reflect a bi-modality in age. 
Hence, changes in the properties of the disk could still be smooth 
as expected from a quiescent internal dynamical evolution such as 
through radial migration of stars. 

Indeed, after deriving the average orbital properties of G dwarf 
stars as function of abundance [a/Fe] vs. metallicity [Fe/H], we 
found that there is a smooth variation with [a/Fe]. At the same time, 
however, the average orbital properties of the a-old stars do change 
with decreasing [Fe/H]: the orbital angular momentum decreases, 
and the orbits become significantly non-circular and thicker. Be- 
cause a basic feature of radial migration is that stars remain on 
near-circular orbits, we next separated the stars based on the circu- 
larity of their orbits, and we found that the Milky Way's thick disk 
is unlikely to have originated through radial migration alone. 

The a-old stars on near-circular orbits show smoothly varying 
properties with [a/Fe], consistent with radial migration into the So- 
lar neighborhood. These stars, which seem to play an important, if 
not dominating role, in the traditional thick disk, dilute any poten- 
tial kinematical or chemical distinction between thin and thick disk. 
Such a smooth extension of the thin disk even brings into question 
whether there is actual a distinct thi ck disk component as implied 
by the photometric decomposition ( iBovvet alJboilbl) . However, 
we also found a-old, mainly metal-poor stars on eccentric orbits 
which are difficult to explain with radial migration, but might have 
formed through early-on gas-rich mergers. Their average orbital 
properties do not show any trends with [a/Fe] nor [Fe/H], and the 
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values are significantly different from those on circular orbits, with 
around twice lower azimuthal rotation velocity, and twice larger 
vertical velocity dispersion. Moreover, as the vertical density pro- 
file of these old stars on eccentric orbits is well described by a sin- 
gle exponential with a scale height of 865 pc, they actually seem to 
form a "traditional" distinct thick disk component. 

To (dis)prove such a mixed origin of the Milky Way's thick(er) 
disk, more direct comparisons with models are essential, which we 
hope to have accommodated by also providing orbital properties 
that can be more easily inferred from the models. At the same time, 
it will be very worthwhile to apply a chemo-orbital analysis as in 
this paper to other types of stars to also extend the observational 
constraints to larger height and radius. 
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